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Abstract 

A superstrong magnetic field stimulates the spontaneous production of positrons 
by naked nuclei by diminishing the value of the critical charge Z cr . The phenomenon 
of screening of the Coulomb potential by a superstrong magnetic field which has been 
discovered recently acts in the opposite direction and prevents the nuclei with Z < 52 
from becoming critical. For Z > 52 for a nucleus to become critical stronger B are 
needed than without taking screening into account. 

1 Introduction 

In a superstrong external magnetic field the Coulomb potential becomes screened [H [2] . The 
screening occurs at the one-loop level; the corresponding Feynman diagram is the polarization 
operator insertion into the photon propagator. This phenomenon leads to the finiteness of 
the ground state energy of a hydrogen atom in the limit of infinite magnetic field B (without 
screening the ground state energy diverges as — In 2 B). In [3] an analytical formula for the 
screened Coulomb potential has been derived. It describes the behavior of the potential 
along the magnetic field and determines the atomic energies. In Sect. II of this paper we 
will derive the formula which describes the behavior of the screened Coulomb potential in the 
direction transverse to the magnetic field. We will see that for B 3> m 2 /e 3 in the transverse 
direction the Coulomb potential is screened at all distances p = y/x 2 + y 2 > 1/ \J e 3 B unlike 
in the longitudinal direction, where the screening takes place at l/m e > z > 1/V e 3 B, and 
in complete analogy with D = 2 QED with light fermions analyzed in 0j . In Sect. Ill we 
will investigate how the effects of higher loops modify the one-loop result for the Coulomb 
potential. The contributions of higher loops are suppressed by powers of the fine structure 
constant a. So with high accuracy the potential is determined by the one-loop result. 

In papers [HE], [3] the spectrum of energies on which the lowest Landau level (LLL) splits 
in the proton electric field was found by solving the corresponding Schrodinger equation. 
According to [3] the ground state energy of hydrogen in the limit of infinite B equals E = 
— 1.7 keV, so the use of the nonrelativistic Schrodinger equation is at least self consistent. 
However, the size an of the electron wave function for B > m 2 /e 3 in the direction transverse 
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to the magnetic field is much smaller than the electron Compton wavelength, a# = 1/ yeB < 
e/m e <C l/m c , which makes the nonrelativistic approach a bit suspicious. That is why 
starting from Sect IV we will study the ground state energy of the electron in a hydrogen-like 
ion in the presence of an external magnetic field by analyzing the Dirac equation. Without 
taking screening into account this problem was considered in paper [5] (see also j6]), soon after 
it was found that a hydrogen-like ion becomes critical at Z = 170: the electron ground level 
sinks into the lower continuum (eq < —m e ) and the vacuum becomes unstable by spontaneous 
e + e~ pairs production. These results were obtained by solving the Dirac equation for an 
electron moving in the field of a nucleus of finite radius. That the phenomenon of criticality 
appears only in the framework of the Dirac equation is an additional motivation to go from 
Schrodinger to Dirac. 

According to [3 |6] the external magnetic field diminishes the value of the atomic charge 
Z cr at which the electron ground level enters the lower continuum. It happens because a 
large magnetic field makes the electron motion quasi-one-dimensional, and in d — 1 the 
potential l/\z\ is more singular than in d = 3 the potential 1/r. 

In Sect. IV from the numerical solution of the Dirac equation for the ground electron 
level of a hydrogen atom in the Coulomb potential we will find that the corrections to the 
nonrelativistic results are small and that the estimate SE = \E® — Ey ch -\ ~ (i?Q Ch ) 2 /m e works 
well. 

In Sect. V we will study how screening modifies the results of paper [5j^J The value of 
the magnetic field B^ at which an ion with charge Z becomes critical increases because of 
screening and only ions with Z > 52 can become critical. 

Let us point out at a major difference between our results and those of [5]. In [5] the 
lower limit on Z above which the ion becomes critical originates from the non-zero size of 
its radius, which "cuts off" the singularity of the Coulomb potential. At the opposite, our 
limit Z > 50 is universal and holds true even for a point-like nucleus. That the ion has a 
finite size only slightly strengthens the constraint. 

Our results are summarized in Table [5] and Fig. [2j in which the value of B/Bq at 
which the ion becomes critical is plotted as a function of its electric charge Z cr . The results 
including screening are given by the blue curve. If one omits screening, one gets the dashed- 
green curve: at the same value of Z cr , criticality with screening is seen to occur at larger 
B, enormously larger for Z = 50. The underlying mechanism that prevents its occurrence 
is the freezing of the ground state energy of the electron, which cannot go into the lower 
continuum. For Z > 20 the screening of the ground state energies occurs in a relativistic 
regime. 

2 3D pattern of the screened Coulomb potential 

Our starting point is the expression for the electric potential in the momentum representation 
<$>(k) where the one-loop contribution to the photon polarization operator in the external 
magnetic field is taken into account. It greatly simplifies when the external magnetic field is 
larger than the Schwinger field: B > B = m 2 e /e (we use Gauss units, e 2 = a = 1/137.03...). 
For such a strong magnetic field the polarization operator is dominated by the contribution 
of electrons which occupy LLL, and a simple and rather accurate interpolation formula for 

1 The necessity of such consideration was stressed in [2]. 
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it was suggested in [4]. With the help of (8) and (12) from [I] we obtain: 



Aire 



(1) 



and 



${z,p) = 47re y 



e fk ^ +ik w z dk ll d 2 k ± /(27r) 



(2) 



where we assume that the magnetic field B is directed along the z axis and k\\ is the mo- 
mentum component parallel to B, while k± and p are vectors in the plane transverse to the 
magnetic field. 

The expression for the electric potential of a pointlike charge in the direction of the 
magnetic field at p = was obtained in 



e 

\z\ 



l-e 



?>™l\z\ + e -^/(2/-K)e^B+&ml\z\ 
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For B <C 3iiml/e 3 the potential equals Coulomb up to small, power suppressed, terms, 
while for B ^> 3itm 2 e /e 3 we get: 
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where Z 



ln( 



The behavior of the potential in the transverse plane (z = 0) can also be found analyti- 



cally in the limit B 3> 3nml/e 3 . Performing the integration in ([2| (for p > an the exponent 
in the denominator can be neglected) we obtain: 



$(0,p) 



exp(- v / (2/7r)e 3 5p) , p < l 



■7rmg e 
e 3 B p 



k < P , 



(5) 



and the Coulomb potential is screened at large p in complete analogy with the D = 2 case, 
see H], Eq. (10). 

' 1 "qi 



For |z| ^> l/m e the values <C m e dominate in the integral (|2j) and we get: 

e 



z»l/m e 



2 2 + (1 + 



e3B 
37rm.2 



)P 2 



(6) 



In Fig. [T] the equipotential lines are shown. The behavior of the screened Coulomb 
potential in the transverse plane was found numerically in [2J, Fig. 2. $(0,p) given by 
expression ^ is close to the results shown in Fig. 2 in [2] for b = B/B = 10 4 -j- 10 6 , within 
the domain of validity of ([5]), and deviate from the curve for b = 10 3 since it corresponds to 
B < 37rmg/e 3 , where ^ is inapplicable. 
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Figure 1: The equipotential lines at B/B Q = 10 4 . The dashed line corresponds to V z 2 + p 2 = 

% /(2/7r)e 3 S V 3irm1 ' 

Finally for 37rm 2 /e 3 ^> B ^> m 2 /e expanding ^ we get (0 is the angle between the two 
vectors f and B): 

${r) = ^i\-^(BlB Q )wo?o\ , r = ^p 2 + z 2 ^> l/m e , (7) 

which coincides with the result obtained in [7] where the expression for the photon polariza- 
tion operator at B > Bq was obtained as well. 

3 Higher loops 

The expression for the screened Coulomb potential was obtained in p]-[3] from the one-loop 
contribution to the photon polarization operator in an external magnetic field B S> Bq. In 
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momentum space it looks like (see (13) in [3]): 



Attc 

k = k ± = 0) = ^- . (8) 

j 2 I 2e3£ fc || 
'|| 7r fe| +6m,2 

If n-loop diagrams contain terms ~ e 3 .B(e 3 -B/A; 2 ) n_1 they would drastically change the 
shape of the potential in coordinate space. 

To calculate the radiative corrections one should use the electron propagator G(k) in an 
external homogeneous magnetic field B. Its spectral representation is a sum over Landau 
levels and for B B the contribution of the lowest level dominates [8] : 

GW = c -^ (1 _ !7l72) |i±^ , (9) 

^0,3 " l e 

where the magnetic field is directed along the third (or z) axis (B = (0,0, B)), fco,3 = 
^o7o ~~ ^373) k± is the component of the momentum normal to the magnetic field and the 
projector (1 — 27172) selects the virtual electron state with spin opposite to the direction 
of the magnetic field. The contributions of the excited Landau levels to G yield a term in 
the denominator proportional to eB and they produce a correction of order e 2 = a in the 
denominator of (|8]). 

Two kind of terms contribute to the polarization operator at the two-loop level. First, 
there are terms in the electron propagators which represent the contributions of higher 
Landau levels. Just like in the one-loop case they produce corrections suppressed by e 2 
in the denominator of i.e. terms of the order e 5 B which can be safely neglected in 
comparison with the leading ~ e 3 B term. Second, there is the contribution from the leading 
term ^ of the electron propagator. Let us consider the simplest diagram: the photon 
dressing of the electron propagator. Neglecting the electron mass we get: 

7 M (1 - 27172)^0,37/, = -2[feo,3 - *&o,3727i] = -2^o,3(l + 27172) , (10) 

which gives zero when multiplied by the external propagator ^ of the electron, since (1 + 
27172) (1 — 27172) = 0. This result is a manifestation of the following well-known fact: in 
D = 2 massless QED (Schwinger model) all loop diagrams are zero except the one-loop 
term in the photon polarization operator. That is why the two loop diagrams in which 
the propagators of virtual electrons are given by ^ give contributions to the polarization 
operator proportional to m 2 for /cjj 3> m 2 while in the opposite limit fcjj m 2 they should 
be proportional to A; 2 in order for the photon to stay massless. So these terms are of the 

order of a(e 3 B) ^m 2 /c 2 /(fc 2 + m 2 ) 2 ) and they are not important. 

The generalization of the above arguments to higher loops is straightforward. 

To conclude this section let us note that an analogous statement about the unimportance 
of the two-loop terms was made in [9]. 
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4 Dirac equation with a screened Coulomb potential, 
Z = 1 



The ground state electron energy of a hydrogen-like ion with electric charge Z in an external 
magnetic field was analyzed in [5] in the framework of the Dirac equation]^] In a strong 
magnetic field (a# = 1/y/eB <C l/(mZe 2 )) the electron spectrum consists of the Landau 
levels splitted into Coulomb sublevels; the ground level belongs to LLL. In complete analogy 
with the nonrelativistic problem the adiabatic approximation is applicable. Averaging over 
the fast motion of the electron in the plane transverse to the magnetic field, the Dirac 
equation for the electron on LLL was reduced in [5] to two first order one-dimensional 
differential equations: 

g z - (e + m e - V)f = 0, ( , 

f z +(e-m e -V)g = 0, 1 ) 

where e is the energy eigenvalue of the Dirac equation; g z = dg/dz, f z = df/dz; the 



bispinor ip e = (^ e ) of the electron is decomposed into 



■r'< ( n(z)e X p(-p 2 /4a 2 H ) ) \' 







i/(z)cxp(-p 2 /4a|-) 

They describe the electron motion in the effective potential V(z): 



V{z) 



-r J nVP 2 + ^ 2 )exp l-£r)pdp , (12) 



where V(r) = —Ze 2 /r, r 2 = p 2 + z 2 . At large distances \z\ 3> an the effective potential 



equals Coulomb, and the solutions of the equations (11) exponentially decreasing at \z\ — > oo 
are linear combinations of Whittaker functions. At short distances the equations (11) can 
be easily integrated for |V"(,z)| 3> \s ± m e \, which is equivalent to the following inequality: 
z Ze 2 /(2m e ). Matching short and large distance solutions at 

Ze 2 /(2m e ) > z > a H (13) 

gives an algebraic equation for the ground state energy (it coincides with Eq. (22) in [5] in 
the limit R/clh <C 1, where R is the nucleus radius): 



Z< In i 2 * ml " \ + arctan { + argT ( f£l + iZe 2 



eB \Vm e -e \ J m 2 _ F 2 



-ar g r(l + 2iZe 2 ) - ^-(ln2 + 7 ) = l + nvr , (14) 

where 7 = 0.5772... is the Euler constant, and the argument of the gamma function is given 
by 

arg T(x + iy) = - jy + Xj (| - arctan x + \ _ 1 ) ' ^ 15 ^ 

For the ground level at e > one should take n = 0, while for e < it should be changed to 
n = —1. 



2 Let us note papers |10j in which the relativistic corrections to the hydrogen spectrum in a strong magnetic 
field are discussed. 



6 



According to (14) when the magnetic field increases the ground state energy goes down 



and reaches the lower continuum. The value of the magnetic field at which this happens is 



determined by (20) (see below) 



A matching point exists only if B 3> <im e /{e(Ze ) ) (see (13)) and (14) is valid only for 



these values of the magnetic field. However, as was checked in [5] from (14) in the non- 
relativistic regime Ze 2 <l,m-£<m, a formula can be deduced which is a valid solution 



of the non-relativistic problem and extends the domain of validity of eq. (14). 



Thus, without taking screening into account, from (14) we can obtain the dependence 



of the ground state energy of a hydrogen atom on the magnetic field for B 3> 4m 2 e /e b . In 
order to find the ground state energy at B < 4m 2 , /e 5 and to take screening into account we 



solve the equations (11) numerically. This system can be transformed into one second order 

1 /2 

differential equation for g(z). By substituting g(z) — (e + m e — V) x( z ) a Schrodinger- 
like equation for the function x( z ) was obtained in [5]r] 



fx 

dz 2 



+ 2m e (E - U)x = 



(16) 



E=*^*, U=^V-^V 2 + V " - + 3/8(y ' )2 



2m e ' m e 2m e 4m e (e + m e — V) m e (e + m e — Vy 



where e is the energy eigenvalue of the Dirac equation and V(z) is given in (12). We in- 



tegrated (16) numerically in the present work. Leaving a detailed discussion for a future 
publication [12J let us only note that, while for z ^> l/m e the last three terms in the expres- 
sion for U are much smaller than the first one (the only one remaining in the nonrelativistic 
approximation), at z < l/m e the relativistic terms dominate and are very big for B 3> Bq 
at z ~ an which makes numerical calculations very complicated. 

In Table [T] the results for the ground state energy of a hydrogen atom without screening 
are presented. The values of the magnetic field in units of B are given in the first column, 
while in columns 2-5 the values of A are given. By definition]^] 



2 2 4 

e — m e __ m e e 2 



E = e - = —A . (17) 

2m P 2 v ; 



From Tabled] we see that: 



1. the results of the numerical integrations of the Schrodinger and Dirac equations coincide 
within four digits: 

• with the analytical Karnakov-Popov formula for the ground state energy (n p = 
m = 0) |T3] E] in the case of the Schrodinger equation; 



with formula (14) for Z = 1 in the case of the Dirac equation; 



2. for the relativistic shift of energy the following estimate works: 



3 This trick was exploited by V.S. Popov for the qualitative analysis of the phenomenon of critical charge. 
Let us note that the definition of A used in [5] differs from our: A0 = e 2 A. 
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Table 1: Values of A for Z = 1 without screening obtained from the Schrodinger and Dirac 
equations. They start to differ substantially at enormous values of the magnetic field. 
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38 


14.88 
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16.93 


16.93 


16.* 


19 


16.89 


10 7 


18.98 


18.98 


18.93 


18.92 


10 8 


21.06 


21.05 


20.98 


20.98 


10 9 


23.16 


23.15 


23.05 


23.05 


10 10 


25.27 


25.27 


25.14 


25.13 


10 11 


27.40 


27.40 


27.23 




10 12 


29.54 


29.54 


29.33 




10 15 


36.03 




35.64 




10 20 


46.99 




46.11 




10 25 


58.07 




56.40 




10 30 


69.22 




66.38 




10 35 


80.43 




75.98 




10 40 


91.67 




85.10 




10 45 


102.95 




93.67 




1Q 50 


114.25 




101.62 




10 55 


125.57 




108. 


89 





E 



Dirac -^Schr 



E\ 



Schr 



Schr 



5A~e 4 A 3 /4 . 

To take screening into account, the following formula for the effective potential should 



be used in (16) instead of (12): 

Ze 2 



V(z) 



I _ e-V^W + e -y/W*)e?B+Gml\z\ J 



°° -p 2 /2a% 



pdp 



(19) 



where Z — 1 for hydrogen. 

The freezing of the ground state energy is due to a weaker singularity of the potential 
with screening (19) at z — > for B — > oo than that of the potential without screening 
(12). While the non-screened potential behaves like 1/z at small z, the screened potential is 
proportional to 5(z) because, when B — > oo, the width of the region where it behaves like 
1/z shrinks to zero [U [2] . 

In Table |2]the results of the analytical formula for A with the account of screening derived 
in [3] for the Schrodinger equation are compared with the results of the numerical integration 
of the Dirac equation. We see that in the case of screening the relativistic shift of energy is 
also very small, and due to it the ground state energies become a little bit higher, just like 
without taking screening into account. The freezing of the ground state energy occurs at 
B/B = 10 3 10 4 , when B « 3nm 2 Je 3 . 



Table 2: Values of A for Z = 1 with screening. 



B/B 


Eq. (57) from [3] 
(Schrodinger) 


Numerical results 
(Schrodinger) 


Numerical results 
(Dirac) 


10° 


5.7 


5.7 


5.7 


10 1 


7.4 


7.4 


7.4 


10 2 


9.1 


9.1 


9.1 


10 3 


10.5 


10.6 


10.6 


10 4 


11.1 


11.2 


11.2 


10 5 


11.2 


11.3 


11.3 


10 6 


11.2 


11.4 


11.3 


10 7 


11.2 


11.4 


11.3 


10 s 


11.2 


11.4 


11.3 



5 Screening versus critical nucleus charge 

According to [5] nuclei with Z > 40 become critical in an external B (for smaller Z the values 
of an at which the criticality is reached become smaller than the nucleus radius, the Coulomb 
potential diminishes and thus the ground level does not reach the lower continuum). 

In Table [3] one can see the dependence of the ground state electron energy Eq on the ex- 



ternal magnetic field for Z = 40. The numerical solutions of (16) are in good correspondence 



with the values of Eq obtained from (14). The numerical results with screening are shown in 
the last column; we see that freezing occurs in the relativistic domain e « — m e /2 and the 
ground level never reaches lower continuum, e > — m e- 

In Table [4] we compare freezing energies for different Z obtained numerically from the 
nonrelativistic Schrodinger equation and from the Dirac equation. We see that for Z > 20 
the freezing occurs in the relativistic regime, where the Schrodinger equation should not be 
used. Let us stress that the value of B at which the freezing occurs does not depend on Z. 



From (14) we obtain in the limiting case E — > — m e an equation which defines the value 



of the magnetic field at which a nucleus with charge Z becomes critical without taking 
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Table 3: Values of Eo/m e for Z = 40. 



B/B 



Eq. (14) 



(Dirac 



Numerical results 
(Dirac) 



Numerical results 
with screening (Dirac) 



10° 
10 1 
10 2 
10 3 
10 4 
2 ■ 10 4 

10 5 
10 6 
10 7 
10 s 



0.819 
0.653 
0.336 
-0.158 
-0.758 
-0.926 
at B/Bq 



2.85 



0.850 
0.667 
0.339 
-0.159 
-0.759 
-0.927 
10 4 , s 



0.850 
0.667 
0.346 
-0.0765 
-0.376 
-0.423 

-0.488 
-0.524 
-0.535 
-0.538 



Table 4: Values of freezing ground state energies for different Z from the Schrodinger and 
the Dirac equations. In order to find the freezing energies we must take B/Bq 3> 37r/e 2 . In 
numerical calculations we took B/B = 10 8 . 



z 


/ t , \ numerical 

(E^\ , keV 

V u J Schr 


/ i \ numerical 

( Er. — m e ) , keV 

V / Dirac 


1 


-1.7 


-1.7 


10 


-88 


-87 


20 


-288 


-273 


30 


-582 


-519 


40 


-966 


-787 


49 




-1003 



screening into account (it coincides with Eq. (32) from [5]): 

B ofy 2 , 2 / , 7r-2argr(l + 2^ cr e 2 ) ^ 

— = 2{Z cr e ) exp -7 + — — . (20) 

-DO \ ^crZ ) 

This equation is used to calculate the numbers in the second column of Table [5j 

From Table [5] we see that with the account of screening only the atoms with Z > 52 
become supercritical at the values of B/B shown in the fourth column. Because of screening 
a larger B is needed for a nucleus to become supercritical and the nuclei with Z < 52 never 
reach supercriticality. This phenomenon is illustrated in Fig. [2] 

From Tables [Tj [3j and [5] we see that ( 14 ) is very good in describing the dependence of 



the energy on the magnetic field; at least a numerical integration produces almost identical 
results. In Table [6] we demonstrate several cases where the accuracy of (fi~4j) is not that 
good. It happens at low B/B since the matching condition B > 4m 2 /(e(Ze 2 ) 2 ) fails and 
when Eq is relativistic. However, B should not be too low to make the adiabaticity condition 
clb S> cih, or B > (Ze 2 ) 2 m 2 /e applicable. 

Textbooks [H] contain detailed consideration of the phenomenon of critical charge. 
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Table 5: Values of B/Bq at which So = —m e according to the Dirac equation and nuclei be- 
come supercritical without (column 2,3) and with (column 4) taking screening into account. 



Z cr 


Eq. (20) 


Numerical results 
without screening 


Numerical results 
with screening 


90 


118 


116 


122 


85 


157 


154 


164 


80 


213 


210 


229 


75 


301 


297 


335 


70 


444 


438 


527 


65 


689 


681 


923 


60 


1144 


1133 


1964 


55 


2068 


2053 


6830 


54 


2357 


2340 


10172 


53 


2699 


2681 


17012 


52 


3107 


3087 


35135 


51 


3594 


3572 


1.20 ■ 10 5 


50 


4181 


4157 


1.14- 10 7 


45 


9826 


9787 




40 


28478 


28408 




35 


1.12 • 10 5 


1.12 ■ 10 5 




30 


6.99 • 10 5 


6.98 • 10 5 




25 


9.27- 10 6 


9.27- 10 6 





Table 6: Values of e /m e at B/B = 5. 



Z 



Eq. (14) 



(Dirac) 



Numerical results 
(Dirac) 



90 
80 
70 
60 
50 
40 
30 
20 
10 
1 



0.2050 
0.3096 
0.4139 
0.5171 
0.6185 
0.7165 
0.8086 
0.8914 
0.9596 
0.998745 



0.2512 
0.3539 
0.4542 
0.5516 
0.6454 
0.7349 
0.8188 
0.8952 
0.9601 
0.998745 



6 Conclusions 

A magnetic field plays a double role in the critical charge phenomenon. By squeezing the 
electron wave function and putting it in the domain of a stronger Coulomb potential it 
diminishes the value of the critical charge substantially [5]. However, for nuclei with Z < 52 
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Figure 2: The values of Bf r : a) without screening according to (20), dashed (green) line; b) 
numerical results with screening, solid (blue) line. The dotted (black) line corresponds to 
the field at which an becomes smaller than the size of the nucleus. 



to become critical such a strong B is needed that the screening of the Coulomb potential 
occurs and acts in the opposite direction: the electron ground state energy freezes and the 
nucleus remains subcritical in spite of growing B. 

We are grateful to O. Babelon, V.L. Chernyak, V.S. Popov, A.I. Rez, V.B. Semikoz, A.E. 
Shabad and S. Teber for useful discussions and remarks. S.G. and M.V. were partially sup- 
ported by the grant RFBR 11-02-00441 and by the grant of Russian Federation Government 
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